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Abstract 

We study the influence of magnetised crusts on torsional shear mode oscillations of magnetars. 

In this context, we employ magnetised crusts whose ground state properties are affected by Landau 

quantisation of electrons. The shear modulus of magnetised crusts is enhanced in strong magnetic 
^ ,- 

fields > 10 G. Though we do not find any appreciable change in frequencies of fundamental 
torsional shear modes, frequencies of first overtones are significantly affected in strong magnetic 
fields. Furthermore, frequencies of torsional shear modes calculated with magnetised crusts are in 
good agreement with frequencies of observed quasi-periodic oscillations. 

PACS numbers: 
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I. INTRODUCTION 



Soft gamma repeaters (SGRs) are characterised by their sporadic and short bursts of soft 
gamma rays. Luminosities in these bursts could reach as high as ~ 10 41 ergs s _1 . There 
are about 5 SGRs known observationally. Evidences of stronger emissions of gamma rays 
from SGRs were observed in several cases. These events are known as giant flares in which 
luminosities are ~ 10 44 — 10 46 ergs s _1 . So far three cases ofgiant flares were reported and 
those are SGR 0526-66, SGR 1900+14 and SGR 1806-20 [JUl. In giant flares, the early 
part of the spectrum was dominated by hard flash of shorter duration followed by a softer 
decaying tail of a few hundreds of seconds. 

SGRs are very good candidates for magnetars which are neutron stars with very high 



surface magnetic fields ~ 10 15 G j^, Giant flares might be caused by the evolving 
magnetic field and its stress on the crust of magnetars. It was argued that starquakes 
associated with giant flares could excite Global Seismic Oscillations (GSOs)[6]. Torsional 
shear modes of magnetars with lower excitation energies would be easily excited. In this 
case, oscillations are restored by the Coulomb forces of crustal ions. Furthermore, the shear 
modes have longer damping times. Quasi-periodic oscillations (QPOs) were found in the 
decaying tail of giant flares from the timing analysis of data (2 - 4 ] . These findings implied 
that QPOs might be torsional shear mode oscillations of magnetar crusts 6|. Frequencies 
of the observed QPOs ranged from 18 Hz to 1800 Hz. 

It was noted from earlier theoretical models of QPOs that the observed frequencies in 
particular higher frequencies could be explained reasonably well using torsional shear os- 
cillations of magnetar crusts (3), [sUlO | . On the other hand, lower frequencies of observed 
QPOs might be connected to Alfven modes of the fluid core. This makes the study of the 
oscillations of magnetar crusts more difficult. There were attempts to explain frequencies 
of QPOs using Alfven oscillations of the fluid core without considering a crust llMl3|. The 
coupling of Alfven oscillations of fluid core with the shear mode oscillations in the solid crust 



due to strong magnetic fields in magnetars was already studied by several groups [14M18 ] . 
It was argued that torsional shear modes of the crust might appear in GSOs and explain 
frequencies of observed QPOs for not very strong magnetic fields despite all these complex 



problems 



19|. 



Nuclear physics of crusts plays an important role on torsional shear modes of magnetar 
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crusts. In particular, the effects of the nuclear symmetry energy on the shear mode fre- 



quencies were investigated recently 



101 ] . It may be worth noting here that torsional shear 



mode frequencies are sensitive to the shear modulus of neutron star crusts. Furthermore, 
the shear modulus strongly depends on the composition of neutron star crusts. Earlier stud- 
ies of torsional shear mode oscillations exploited only non-magnetic crusts. Magnetic fields 
in magnetars may influence the ground state properties of neutron star crusts. Recently, 
we have investigated the influence of Landau quantisation of electrons on the compositions 
and equations of state 



those properties 
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CoS) of outer and inner crusts and obtained appreciable changes in 
21] . This, in turn, might influence the shear modulus of crusts and 
torsional shear mode oscillations in magnetars. This motivates us to study torsional shear 
mode oscillations of magnetars using magnetised crusts. 

We organise the paper in the following way. We describe models for torsional shear 
mode oscillations, shear modulus and compositions and EoS of magnetised crusts in Sec. 
II. Results of this calculation are discussed in Sec. III. Section IV gives the summary and 
conclusions. 



II. FORMALISM 



Earlier calculations of torsional shear mode oscillations were performed in Newtonian 



gravity 



23] as well as in general relativity 



an 
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25] with and without 



magnetic fields. In many of these calculations, the magnetised crust was decoupled from the 
fluid core. 

As we are interested in the effects of magnetised crusts on torsional shear mode frequen- 
cies, we consider a free slip between the crust and core. Here we calculate torsional shear 
mode frequencies following the model of Refs.^, [25)]. In this case, we study torsional shear 
mode oscillations of spherical and non-rotating relativistic stars in the presence of a dipole 
magnetic field. We neglect the deformation in the equilibrium star due to magnetic fields 
~ 10 16 G because its effects would be smaller than the direct influence of magnetic fields 
on shear mode frequencies. Therefore, we assume the magnetised star to be spherically 
symmetric. The metric used to determine equilibrium stellar models has the form, 



ds * = _ e 2Z dt 2 + e 2A dr 2 + r 2 ^2 + 



(1) 
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The equilibrium models are obtained by solving the Tolman-Oppenheimer-Volkoff (TOV) 
equation with a perfect fluid EoS. The stress-energy tensor of a magnetised relativistic star 
in equilibrium is written as 

where the stress-energy tensor of perfect fluid is 



Tfi = (e + P)uV + Pg^ 



and that of the electromagnetic field is 



rpfJLV 
1 EM 



— \B 2 u^u u + -g^B 2 -B^B v ) , 
4tt V 2 y J ' 



(3) 



(4) 



where the four velocity of the fluid element is = (e - *, 0, 0, 0) and u^u^ — — 1. 

Here we consider an axisymmetric poloidal magnetic field generated by four current = 
(0,0,0, J$). After expanding the four-potential into vector spherical harmonics as = 
ae m (r)sin6d0Pe m (cos6) and similarly for J M , it follows from Maxwell's equations for a dipole 
magnetic field i.e. £ m = 1 



_ 2A d a\ ^ ^ yY')e~ 2A 1 



dr 2 



dr 



2a x 



-47TJi 



(5) 



where prime denotes derivative with respect to r. The exterior solution for ai is given by 

3/i fe r 2 



(ex) 



3M 3 
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2M 



2M 2M 2 



r 



(6) 



where //& is the magnetic dipole moment for an observer at infinity and M is the mass of 
the star. The interior solution can be obtained by integrating Eq. ([3]) assuming a current 



distribution ji as was done by Ref. 



9|. The components of the magnetic field are 
2e A cos6 



By 

Bo 



-ax 



~ A sin6 



dax 
dr 



(7) 



The perturbed equations describing shear mode oscillations are obtained by linearising 
the equations of motion of the fluid and the magnetic induction equation 0, . Torsional 
shear modes are incompressible and do not result in any appreciable density perturbation in 
equilibrium stars. Consequently, one may adopt the relativistic Cowling approximation and 



neglect metric perturbations Sg, lu =0 



261 ] . Torsional shear modes are the results of material 



velocity oscillations. These modes are odd parity (or axial type) modes. We consider axial 
type perturbation in the four velocity and the relevant perturbed matter quantity is the 
0-component of the perturbed four velocity du^ ^] 



du* = e-*d t y(t,r)—-d Pi(cos9) 



sinO 



(8) 



where d t and d$ correspond to partial derivatives with respect to time and 9, respectively, 
Pi(cos8) is the Legendre polynomial of order I and y(t, r) is the angular displacement of the 
matter. It is to be noted that the radial and angular variations of azimuthal displacement of 
stellar matter lead to shears of the crystal lattice in neutron star crusts which are described 



by the shear tensor S 



24| . Further, the shear stress tensor is given by T^ u = —2fiS l 



(Wl 



where \x is the isotropic shear modulus. The linearised equations of motion includes the 

contribution of ST^ [9j . 

Assuming a harmonic time dependence for y(t, r) = e 1U3t ytfj) and neglecting I ± 2 terms, 
one obtains the eigenvalue equation for the mode frequency 



M+(l + 2A 1 )^ T 
7rr* 



y" 



{ (~ + *' - A') n + p! + (1 + 2Ai)^ [(*' - A') oi + 2ai'] } V 



e + p + (l + 2Ai) 



2 \ \ /2 

a l \ ,2A A iai 

~ 2nr 2 



(A -2) 



Aiai 



2?rr 4 



(9) 



where A = £(£ + 1) and Ai = -£{£ + l]/(2£ - 1)(2£ + 3). Equation © reduces to the 
non-magnetic case when we put ai = 19|. With suitable choice of new variables, Eq.()9]) 
results in a system of first order ordinary differential equations [9] . As torsional shear modes 
are confined to the crust in our calculation, we impose a zero traction boundary condition 
at the interface between the core and the crust as well as the zero torque condition at the 
surface 0,13]. These conditions imply y = at surface (r = R) of the star and the interface 
(r = R c ) of the crust and core. Finally, we estimate eigenfrequencies by solving two first 
order differential equations. 

The knowledge of the shear modulus of the magnetised crusts is an important input in 
the eigenvalue equation (Eq.Q) for the shear mode calculation. Here we adopt the following 
expression of the shear modulus at zero temperature 271, [28] 

MZef 



fi = 0.1194- 



(10) 
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where a = 3/(47mj), Z is the atomic number of a nucleus and n» is the ion density. This 
form of the shear modulus was obtained by assuming a bcc lattice and performing direc- 
tional averages [29( . Further the dependence of the shear modulus on temperature was also 
investigated with the Monte Carlo sampling technique by Strohmayer et al. [28|]. The com- 
position and equation of state of neutron crusts are essential ingredients for the calculation 
of the shear modulus as it is evident from Eq. fllOl) . In the following paragraphs, we describe 
the ground state properties in outer and inner crusts in presence of strong magnetic fields. 

The outer crust is composed of nuclei immersed in a uniform background of a nonin- 
teracting electron gas. Neutrons start coming out of nuclei when the neutron drip point 
is reached. This is the beginning of the inner crust where nuclei are placed both in free 
neutrons as well as electrons. Nuclei are arranged in a bcc lattice in neutron star crusts. 
The Wigner-Seitz approximation is adopted in this case. Each lattice volume is replaced by 
a spherical cell with one nucleus at its center. The cell contains equal numbers of protons 
and electrons so that it is charge neutral. Moreover, there is no Coulomb interaction among 
cells. We obtain an equilibrium nucleus in the outer crust minimising the Gibbs free energy 
per particle at a fixed pressure varying mass and atomic numbers 30j 

Etot + P 



n b 



(11) 



where rib is the baryon density and the energy density E tot = un(Wn + W£) + e e includes 
contributions from the energy of the nucleus (Wn), lattice energy (Wl) of the cell involving 
the finite size effects and free electron gas (e e ) 20] . Similarly, the total pressure is given by 
sum of the pressure due to the lattice and that of the electron gas. We use experimental 
nuclear masses from the atomic mass table 3JJ whenever it is available. Otherwise, the 
theoretical extrapolation is adopted in this calculation 32]. 

We describe the ground state properties of matter of the inner crust using the Thomas- 
Fermi model. The spherical cell that contains neutrons and protons does not define a nucleus. 
We adopt the procedure of Bonche, Levit and Vautherin to subtract the free neutron gas 



of the cell and obtain the nucleus 



33 



341 ] . Under this assumption, the thermodynamic 



potential (flw) of the nucleus is obtained as [33l. |34| 



f2 /v — Q 



NG 



fir 



(12) 



where Qng is the thermodynamical potential for the nucleus plus the free neutron gas and 
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Qg is that of the free neutron gas. The thermodynamical potential is defined as 



(13) 



q=n,p 



where J 7 , fi q and n q are the free energy density, baryon chemical potential and number 
density, respectively. The free energy is given by 



T{n q , Y p ) = J[H + e c + e e ]dr , 



(14) 



where Y p is the proton fraction. Equation (j!4p includes contributions from the nuclear 
ener gy d ensity (%), Coulomb energy density (e c ) and energy density of free electron gas 
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interaction 



35] . The nuclear energy density is calculated using the Skyrme nucleon-nucleon 
38]. We obtain the density profiles of neutrons and protons for the nucleus 



36438 



plus neutron gas as well as neutron gas phases by minimising thermodynamical potentials 



2l|. 



and hence calculate mass and atomic numbers of nuclei of the inner crust 

Now we focus on neutron star crusts in strongly quantising magnetic fields. In particular, 
we showed earlier that the Landau quantisation of electrons strongly influenced ground state 
properties of neutron star crusts in strong mag netic fields ~ 10 17 G [20|, [2lJ. Energy and 
number densities of electrons are affected by the phase space modifications due to Landau 
quantisation of electrons. We immediately note that the electron energy density in Eqs.l fTTl) 
and (fl~4"|) are modified in presence of magnetic fields 20, |2J| . Further, we have to take into 
account the change in the average electron chemical potential due to magnetic fields in the 
outer and inner crusts. It is to be noted that protons are only influenced by magnetic fields 
through the charge neutrality condition. 



III. RESULTS AND DISCUSSIONS 



We investigated the composition and EoS of ground state matter in neutron star crusts in 



strong magnetic fields 



20. 



21] . We noted that the electron number density in the outer crust 



was enhanced compared with the field free case when a few Landau levels were populated 
for magnetic fields > B = 4.414 x 10 16 G [20]. It was observed that this enhancement grew 
stronger when only the zeroth Landau level was populated at a magnetic field strength of 
4.414 xlO 17 G. Consequently, we found modifications in the sequence of equilibrium nuclei 
which was obtained by minimising the Gibbs free energy per nucleon of Eq. ffTTT) . It was 



noted that some new nuclei such as fgSr and 4g 8 Pd appeared and some nuclei such as ||Ni 



and ggNi disappeared in a magnetic field of B = 4.414 x 10 16 G 20] when we compared this 
with the zero field case. The maximum density upto which an equilibrium nucleus would 
exist, increased as the strength of magnetic field increased. This implied that the neutron 
drip point was shifted to higher density in presence of a strong magnetic field with respect 

2o|. We perf 



to the field free case 



nucleon-nucleon interaction 
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brmed the calculation of the inner crust using the SkM 
351 ] . In this case too, we calculated the equilibrium nucleus 
at each density point. Like the outer crust in strong magnetic fields, the electron number 
density was enhanced due to the electron population in the zero Landau level for magnetic 
fields > 10 17 G which , in turn, led to a large proton fraction because of charge neutrality. 
For magnetic fields > 10 17 G, equilibrium nuclei with larger mass and atomic numbers were 



found to exist in the crust 



2l|. Furthermore, the free energy per nucleon of the nuclear 



system was reduced in magnetic fields compared with the corresponding case without a 
magnetic field. 

We calculate the shear modulus using Eq. tflOl and the above mentioned models of mag- 
netised crusts. For this purpose, we have to know the profiles of pressure, energy density 
and shear modulus in a neutron star. Those profiles are obtained by solving the TOV equa- 
tion. In this context, we construct an EoS of dense nuclear matter in strong magnetic fields 
in neutron star core using a relativistic mean field model with the GM1 parameter set as 
described in Ref. 39l44l]. This EoS of dense nuclear matter is matched with the EoS of the 
crust and used in the TOV equation. Figure 1 displays the shear modulus as a function of 
mass density for a neutron star of 1.4 M Q . The shear modulus (dahsed line) corresponding 
to B* = B/B c = 10 3 where B c = 4.414 x 10 13 G, does not show any appreciable change from 
that of the zero field (dotted line) because of large numbers of Landau levels are populated 
in this field. As the field strength is increased, less numbers of Landau levels are populated. 
For B* = 10 4 i.e. 4.414 x 10 17 G, the shear modulus is enhanced due to the population 
of all electrons in the zeroth Landau level. In all three cases, the shear modulus increases 
with mass density well before the crust-core interface. The shear modulus and shear speed 
v s = (/i/p) 1 / 2 are extrapolated to the zero value at the crust-core interface for magnetised as 
well as non-magnetised crusts. It was noted that the shear modulus was found to decrease 
smoothly to zero with increasing density when the pasta phase was considered in the crust 
in absence of magnetic fields 19] . We generate a profile of the shear modulus as a function 
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of radial distance in a neutron star for calculating frequencies of torsional shear modes. 

Now we study the dependence of torsional shear mode frequencies on the compositions and 
the shear modulus of magnetised crusts. Earlier all calculations were performed using non- 
magnetic neutron star crusts. Here we exploit models of non-magnetic as well as magnetic 
crusts which are already described in this section. We consider torsional shear modes of 
a neutron star with 1.4 M & . Frequencies of fundamental (n = 0, £ = 2) torsional shear 
modes are plotted with magnetic fields in Fig. 2. It is observed that the frequency increases 
with magnetic field. However, there is no differences between our results with and without 
magnetised crusts. Figure 3 shows frequencies of torsional shear modes corresponding to 
n = plotted as a function of I values for a 1.4 M neutron star and magnetic field = 10 3 
i.e. 4.414 x 10 16 G. Again we calculate frequencies using the non-magnetic and magnetic 
crusts in this case. The frequency increases with higher I values, but we can not distinguish 
between the results of non-magnetic and magnetic crusts. 

We continue our investigation on frequencies of first overtones (n = 1) of torsional shear 
modes in presence of a magnetic field. Frequencies of first overtones are plotted with i values 
for a neutron star of 1.4 M Q and magnetic field B* = 10 3 i.e. 4.414 x 10 16 G in Fig. 4. 
In this case, the results (dashed line) obtained with magnetised crusts deviate from those 
of non-magnetised crusts (solid line). We obtain qualitatively similar results for B* = 10 4 . 
It was noted that the ratio of the crust thickness to the radius of a neutron star could be 
obtained from the frequencies of overtones [9|. Results of Fig. 4 point to the fact that the 
ratio of the crust thickness to the radius might be influenced by strong magnetic fields. 

The dependence of frequencies of the fundamental mode and higher harmonics on neutron 
star masses is demonstrated in Fig. 5. Here the frequencies corresponding to n = and 
I = 2, 3, 4 are shown as a function of neutron star masses for a magnetic field B = 8 x 10 14 G. 
For each case, frequencies of torsional shear modes decrease with increasing mass, whereas 
higher i values lead to higher frequencies. Frequencies of observed QPOs might put a strong 
constraint on the EoS if masses of neutron stars are known accurately. 

Finally, we compare the calculated frequencies of torsional shear modes with frequencies 
of observed QPOs. This comparison is shown in Table 1. Our results are obtained using 
the magnetised crusts of a 1.2 M & neutron star and estimated magnetic fields B = 8 x 10 14 
G corresponding to SGR 1806-20 and 1.4 M neutron star and B = 4 x 10 14 G for SGR 
1900+14 [42H44j . For SGR 1806-20, calculated frequencies 93 Hz corresponding to n = 
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and t = 12 and above are in very good agreement with the observed frequencies. However, 
lower frequencies 18, 26 and 29 Hz can not be explained with our calculation because those 
are so close that they can not be matched using different harmonics of n = mode. On the 
other hand, the overall agreement of calculated frequencies with the observed frequencies of 
SGR 1900+14 is excellent. 



IV. SUMMARY AND CONCLUSIONS 



We have estimated frequencies of torsional shear modes of magnetars assuming a dipole 
magnetic field configuration. Frequencies are computed using our models of magnetised 
crusts. The shear modulus of magnetised crusts is found to be enhanced in strong magnetic 
fields ~ 4.414 x 10 17 G because electrons populate the zeroth Landau level. It is observed 
that frequencies of fundamental (n = 0, £ = 2) torsional shear modes are not sensitive to this 
enhancement in the shear modulus in strong magnetic fields. On the other hand, frequencies 
of first overtones (n = 1) of torsional shear modes in presence of strongly quantising magnetic 
fields are distinctly different from those of the field free case. Consequently, this might impact 
the ratio of the crust thickness to the radius of a magnetar. We have compared our results 
with frequencies of observed QPOs and found good agreement. Observed frequencies could 
constrain the EoS of magnetised neutron star crusts if masses of neutron stars are known. 
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TABLE I: Frequencies of torsional shear modes are recorded here. Frequencies obtained in the 
calculation using the magnetised crusts 



20, 



211 ] are compared with observed frequencies in SGR 
1806-20 and SGR 1900+14. Our results are calculated using a 1.2 Mq neutron star and magnetic 
field 8 x 10 14 G for SGR 1806-20 and a 1.4 M© neutron star and magnetic field 4 x 10 14 G for SGR 
1900+14. 
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FIG. 1. Shear modulus is plotted as a function of mass density for a neutron star of 1.4 M@ and 
different magnetic fields. 
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Fig. 2. Frequency of fundamental (n = 0, I = 2) torsional shear mode for a neutron star of 1.4 
Mq is shown as a function of magnetic field B* = B/B c where B c = 4.414 x 10 13 G. Results of our 
calculations with and without magnetised crusts are shown here. 
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Fig. 4 Frequency of first overtone (n = 1) of torsional shear mode of £ values with and without 
magnetic crusts of a 1.4 M & neutron star for magnetic field = 10 3 . 
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FIG. 5. Frequency of torsional shear mode corresponding to n = and I = 2, 3, 4 as a function of 
neutron star mass for a magnetic field B = 8 x 10 14 G. 
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